Phase Transition in a Conserved-Mass Model of Aggregation and Dissociation 
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We introduce a new model of aggregation of particles where in addition to diffusion and aggrega- 
tion upon contact, a single unit of mass can dissociate from a conglomerate. This dissociation move 
conserves the total mass and leads to a striking behaviour in the steady state. As the parameters 
are varied, the system undergoes a dynamical phase transition in all dimensions. In one phase the 
mass distribution decays exponentially for large mass whereas in the other phase there is a power 
law distribution of masses which coexists with an infinite mass aggregate. The model is investigated 
analytically within mean field theory, and numerically in one dimension. 



Conservation laws satisfied by the dynamics are known 
to modify drastically the time-dependent behaviour 
of systems which are in thermal equilibrium. There is, 
however, a wide variety of inherently nonequilibrium sys- 
tems in nature whose steady states are not described by 
the Gibbs distribution, but are determined by the dy- 
namics itself. Examples include systems exhibiting self- 
organized criticality several reaction-diffusion sys- 
tems 1^ and fluctuating interfaces Q . Then the question 
naturally arises: What is the role played by conservation 
laws in selecting the steady state of such a nonequilib- 
rium system? The question is particularly significant be- 
cause of the propensity of many nonequilibrium systems 
to organize themselves into critical states, which are espe- 
cially sensitive to such conditions. How do conservation 
laws modify the power laws characteristic of such criti- 
cal states? Can they bring in completely new types of 
behaviour? 

We address these questions for the steady states of 
an important class of nonequilibrium processes, namely 
those involving the twin phenomena of aggregation and 
diffusion. These processes are ubiquitous in nature, and 
arise in a variety of physical settings, for instance, in the 
formation of colloidal suspensions ^ and polymer gels H] 
on the one hand, and aerosols and clouds on the other. 
A recent interesting and important result in this area due 
to Takayasu and collaborators which has widespread 
applications, is that constant injection of single particles 
into such a system leads to a power law distribution of 
particle masses in the steady state. This injection of par- 
ticles from outside of course violates mass conservation. 
In this Letter, we show that a conserved-mass system 
in which injection is replaced by dissociation, exhibits 
strikingly different and even more interesting behaviour: 
the new dissociation moves that conserve the total mass 
induce a novel dynamical phase transition. As the pa- 
rameters of the system are changed, there is a transition 
from one steady state where the mass distribution decays 
exponentially to another where it decays as a power law, 
and in addition develops an infinitely large aggregate. 

The conserved-mass aggregation model (CMAM) dis- 
cussed here has connections to models of gelation Q and 
to the Takayasu model (TM) of particle injection alluded 
to above pi . However it turns out to have a steady state 



structure which is richer than either. In experiments 
as well as theoretical models of irreversible gel forma- 
tion, the steady-state mass distribution P[m) is trivial, 
as there is only a single infinite aggregate, though the 
kinetics of the approach to this state is quite interest- 
ing H]. In the TM on the other hand the injection move 
completely destroys the infinite aggregate and the steady 
state mass distribution P{m) decays as a power law to^"^ 
for large mass, where the exponent r depends on the 
spatial dimension. 

Our conserved-mass model differs from conventional 
models of aggregation in that single particles are allowed 
to chip off from more massive conglomerates. This move 
corresponds to the physical process of single functional 
units breaking off from larger clusters in the polymeriza- 
tion problem. It leads to a replenishment of the lower end 
of the mass spectrum, and competes with the tendency of 
the coalescence process to produce more massive aggre- 
gates. The result of this competition is that two types 
of steady states are possible, and there is a dynamical 
phase transition between the two. In one state, P{m) 
decays exponentially, while the other is more interesting: 
P(m) decays as a power law for large m but in addition 
develops a delta function peak at m = oo. Physically this 
means that an infinite aggregate forms that subsumes a 
finite fraction of the total mass, and coexists with smaller 
finite clusters whose mass distribution has a power law 
tail. In the language of sol-gel transition, the infinite ag- 
gregate is like the gel while the smaller clusters form the 
sol. However, as opposed to the models of irreversible 
gelation where the sol disappears in the steady state, 
in our model the sol coexists with the gel even in the 
steady state. Interestingly, the mechanism of the forma- 
tion of the infinite aggregate in the steady state resem- 
bles Bose-Einstein condensation (BEC), though the con- 
densate (the infinite aggregate here) forms in real space 
rather than momentum space as in conventional BEC. 

Our model can be considered as the conserved counter- 
part of the non-conserved TM The injection move in 
TM that violates the mass conservation is replaced in our 
model by the dissociation move that conserves the mass. 
Besides being a simple model having a self-organized crit- 
ical state, TM has found widespread applications includ- 
ing modeling of river networks and stress distribution 
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in granular media pl| . It is therefore not unreasonable to 
expect that the conserved mass model discussed here will 
also find applications in a wide variety of physical pro- 
cesses. As an example, it can be considered as a simple 
model of river networks in a basin where there is negligi- 
ble rainfall (injection) but small rivulets can break off a 
stream (dissociation). 

The CMAM is defined as follows. For simplicity we de- 
fine the model on a one dimensional lattice with periodic 
boundary conditions although generalizations to higher 
dimensions are quite straightforward. Beginning with a 
state in which the masses are placed randomly, a site is 
chosen at random. If it contains one or more than one 
particle, then one of the following events can occur: 

1. Diffusion and Aggregation: With probability pi, 
the mass at site i moves either to site i — 1 or 
to site i + If it moves to a site which already has 
some particles, then the total mass just adds up. 

2. Chipping (single-particle dissociation) : With prob- 
ability p2, a bit of the mass at the site "chips" off, 
i.e. a single particle leaves site i and moves with 
equal probability to one of the neighbouring sites 
i — I and i + I. 

3. With a probability l~pi~p2, the site is left undis- 
turbed. 

If the site chosen is empty, then nothing happens. The 
same rules hold even if we choose a site with only a sin- 
gle particle, which means that the probability for a single 
particle to move left or right in this model is pi +p2 ■ Note 
the difference with TM: in TM, the move 2 is replaced 
by addition of unit mass to every site with probability 1. 

We first analyze the model within the mean field ap- 
proximation, ignoring correlations in the occupancy of 
adjacent sites. Then we can directly write down equa- 
tions for P(to, t), the probability that any site has a mass 
771 at time t. 



dP{m,t) 
di 

dP{0,t) 
Jt 



= -(pi +P2)[1 + q{t)]P{m,t) +p2P{m+l,t) 
+ P2q{t)P{ni-l,t)+piP*P; m > 1 (1) 

-(pi + P2)q{t)P{0, t) + P2P{l,t) + p^qit). (2) 



Here q{t) = 1 — P(0, t) is the probability that a site is oc- 
cupied by a mass and P*P = X)m'=i P{''n', t)P{m—m', t) 
is a convolution term that describes the coalescence of 
two masses. 

The above equations enumerate all possible ways in 
which the mass at a site might change. The first term in 
Eq. (|l|) is the "loss" term that accounts for the probabil- 
ity that a mass m might move as a whole or chip off to 
either of the neighbouring sites, or a mass from the neigh- 
bouring site might move or chip off to the site in consider- 
ation. The probability of occupation of the neighbouring 



site, q{t) = J2m=i multiphes P{m,t) within the 

mean-field approximation where one neglects the spatial 
correlations in the occupation probabilities of neighbour- 
ing sites. The remaining three terms in Eq. (^ are the 
"gain" terms enumerating the number of ways that a site 
with mass ml ^ m can gain the deficit mass m^m' . The 
second equation Eq. (|2|) is a similar enumeration of the 
possibilities for loss and gain of empty sites. Evidently, 
the mean field equations conserve the total mass. 

To solve the equations, we compute the generating 
function, Q{z,t) = ^("^' *)-^™ from. Eq. (|l|) and 

set dQ/dt — in the steady state. We also need to use 
Eq. @) to write P(l,t) in terms of q{t). This gives us a 
quadratic equation for Q in the steady state. Choosing 
the root that corresponds to Q{z = 0) = 0, we find 



Q{z) 
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where w — P2/P1 and zi^2 = (w + 2^ 2\/w -|- l)/wq. The 
value of the occupation probability q is fixed by mass 
conservation which implies that ^777,^(777) — M/L = p. 
Putting dzQ{z = 1) = p, the resulting relation between 
p and q is 



2p = w{l - q) - wqy/ (zi - l)(z2 - 1) 



(4) 



The steady state probability distribution P{ra) is the 
coefficient of in Q{z) and can be obtained from Q{z) 
in Eq. (||) by evaluating the integral 
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over the contour Co encircling the origin. The singular- 
ities of the integrand govern the asymptotic behaviour 
of P(m) for large 777. Clearly the integrand has branch 
cuts at z = zi^2- For fixed w, if one increases the density 
p, the occupation probability q also increases as evident 
from Eq. (^. As a result, both the roots zi^2 start de- 
creasing. As long as the lower root zi is greater than 1, 
Eq. ^ is well defined and the analysis of the contour 
integration around the branch cut z = zi , yields for large 



P(777) 
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where the characteristic mass, 777* = l/log(zi) and di- 
verges as ~ {qc — q)~^ as g approaches qc = (w -|- 2 — 
2\/w -\- l)/u'. qc is the critical value of q at which zi = 1. 
This exponentially decaying mass distribution is the sig- 
nature of "disordered" phase which occurs for q < qc or 
equivalently from Eq. (0) for p < pc = ^/'w + 1 — 1. 
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When p = pc, we have zi = 1, and analysis of the 
contour around z — zi — 1 yields a power law decay of 
P{m), 



P{m) 



-5/2 



(7) 



As p is increased further beyond pc, q can not increase 
any more because if it does so, the root zi would be less 
than 1 (while the other root Z2 is still bigger than 1) and 
Eq. would be undefined. The only possibility is that q 
sticks to its critical value qc or equivalently the lower root 
zi sticks to 1. Physically this implies that adding more 
particles does not change the occupation probability of 
sites. This can happen only if all the additional particles 
(as p is increased) aggregate on a vanishing fraction of 
sites, thus not contributing to the occupation of the oth- 
ers. Hence in this "infinite-aggregate" phase P{m) has 
an infinite-mass aggregate, in addition to the power law 
decay m~^/^. Concomitantly Eq. (^) ceases to hold, and 
the relation now becomes 



P 



w , 
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where poo is the fraction of the mass in the infinite aggre- 
gate. The mechanism of the formation of the aggregate is 
reminiscent of Bose Einstein condensation. In that case, 
for temperatures in which a macroscopic condensate ex- 
ists, particles added to the system do not contribute to 
the occupation of the excited states; they only add to the 
condensate, as they do to the infinite aggregate here. 

Thus the mean field phase diagram (Fig. 1) of the sys- 
tem consists of two phases, "disordered" and "infinite- 
aggregate" , which are separated by the phase boundary, 
Pc = \/w + 1 — 1. While this phase diagram remains 
qualitatively the same even in 1-d, the exponents char- 
acterizing the power laws are different from their mean 
field values (see Fig. 1). 

We have studied this model using Monte Carlo simula- 
tions on a one-dimensional lattice. Although we present 
results here for a relatively small size lattice, L = 1024, 
we have checked our results for larger sizes as well. 
We confirmed that all the qualitative predictions of the 
mean-field theory remain true, by calculating P(m) nu- 
merically in the steady state. Figure 2 displays two nu- 
merically obtained plots of P{m). The existence of both 
the disordered (denoted by +) and the infinite-aggregate 
phase (denoted by x) is confirmed. In particular, the sec- 
ond curve shows clear evidence of a power-law behaviour 
of the distribution, which is cut off by finite-size effects, 
and for an 'infinite' aggregate beyond. We confirmed that 
the mass Magg in this aggregate grows linearly with the 
size, and that the spread SMagg grows sublinearly, im- 
plying that the ratio 5Magg/Magg approaches zero in the 
thermodynamic limit. The exponent r which character- 
izes the finite-mass fragment power law decay is numeri- 
cally found to be 2.33 ± .02. The difference of this value 



from the mean-field value 2.5 is presumably due to the 
neglect of spatial correlations within mean-field theory. 

We note that in TM, the exponent ttm that charac- 
terizes the power law decay of mass distribution in the 
steady state has an exact value 4/3 in 1-d and 3/2 within 
mean field theory . In the CMAM, we find that in the 
aggregate phase tcmam ~ 2.33 in 1-d and 5/2 within the 
mean field theory. It is therefore tempting to conjecture 
that TCMAM = TTM + although we have no proof of 
this. 

Mass conservation evidently affects the steady state in 
this class of nonequilibrium models, but what are the 
other factors which determine the universality classes? 
We addressed this question by studying the effect of di- 
rectionality on the motion of the masses, leading to a 
finite mass current. For the mass-nonconserved TM, it 
is known that making such a change has no effect on the 
scaling properties Our numerical study of the di- 
rected mass-conserved case shows that directionality 
in fact changes the universality class. As in the undi- 
rected case, the model continues to show a phase tran- 
sition between two phases, without and with an infinite 
aggregate, respectively. However, the exponent r char- 
acterizing the power law decay in the infinite-aggregate 
phase, is different in this case (Fig. 3). In this model, 
r ~ 2.05 within numerical error [|l3|. Clearly, P{m) 
should decay faster than m^^ for large m to keep the 
total mass ^mP(m) finite. Perhaps P{m) decays as 
with additional logarithmic factors. 

Interestingly, the CMAM can be mapped onto a driven 
diffusive lattice gas model, and thence to a model of in- 
terface dynamics, so that our results for aggregation phe- 
nomena have wider significance. A configuration of the 
CMAM is mapped onto a particle-hole configuration as 
follows: Every site with m particles in the aggregation 
model is mapped to a cluster of m successively occupied 
sites, with a vacancy at the rightmost edge in the parti- 
cle model. An empty site maps to an empty site. Thus 
a lattice of L sites in the aggregation-dissociation model 
maps to one of L -I- Np sites in the hard-core particle 
model, where Np is the total number of particles. The 
dynamics defined in the earlier section now translates to 
the following rules for particle motion. If a randomly 
selected site is occupied, then with probability pi, the 
entire particle cluster to which the occupied site belongs 
is moved one site to the left or right. With probability 
P2, the rightmost or leftmost particle from the cluster 
dissociates from the cluster and moves one site to the 
right or left. With probability 1 — pi — P2, the system is 
left undisturbed. The mapping to an interface model fol- 
lows a standard route 1 12|: a particle is mapped to a unit 
segment with a positive slope, while a hole maps to a seg- 
ment of unit length with a negative slope. The equivalent 
interface then allows both for single corner flips at hills 
and valleys, and also for slice- wise moves of connected 
segments of up and down slopes. Under this dynam- 
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ics, interface profiles can develop very different spatial 
structures as compared to those described in customary 
growth processes In the exponential phase, parti- 
cle cluster lengths decay exponentially as is usual, but 
in the infinite-aggregate phase, a single macroscopically 
long stretch with positive slope develops. This unusual 
behaviour differs from descriptions of mound formation 
, in that the infinite stretch is accompanied by the 
formation of fairly long finite stretches with a power-law 
distribution of stretch lengths. This defines a new uni- 
versality class in interface dynamics. 

In conclusion, we point out that several questions still 
remain open. Amongst these is the understanding of spa- 
tial and temporal correlations, and the role of direction- 
ality and dimensionality in influencing the scaling be- 
haviour, in particular, the determination of the upper 
critical dimension. 

We thank Deepak Dhar for useful discussions. 
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Figure Captions 

Fig. 1. The phase diagram in the p-w plane. The dashed line shows the mean field phase boundary. The points 
denoted by (•) are obtained numerically in one dimension. 

Fig. 2. The mass distribution P{m) vs. m for the undirected model in 1-d on a log-log plot for w = 1.0 and p = 0.2 
(shown by x symbols) and for w = 1.0 and p = 3.0 (shown by -f symbols). 

Fig. 3. The mass distribution P{m) vs. m (only the power law part) on a log-log plot for both undirected (+) and 
directed (x) case in the "Infinite-Aggregate" phase at density p = 12.0 and w = 1.0. The straight line with slope —2 
is drawn as a guidance to the eyes. 
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